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Abstract. - Finescale probes of many complex electronic materials have revealed a non- 
uniform nanoworld of sign-varying textures in strain, charge and magnetization, forming me- 
andering ribbons, stripe segments or droplets. We introduce and simulate a Ginzburg-Landau 
model for a structural transition, with strains coupling to charge and magnetization. Charge 
doping acts as a local stress that deforms surrounding unit cells without generating defects. 
This seemingly innocuous constraint of elastic 'compatibility', in fact induces crucial anisotropic 
long-range forces of unit-cell discrete symmetry, that interweave opposite-sign competing strains 
to produce polaronic elasto- magnetic textures in the composite variables. Simulations with ran- 
dom local doping below the solid-solid transformation temperature reveal rich multiscale textur- 
ing from induced elastic fields: nanoscale phase separation, mesoscale intrinsic inhomogeneities, 
textural cross-coupling to external stress and magnetic field, and temperature-dependent per- 
colation. We describe how this composite textured polaron concept can be valuable for doped 
manganites, cuprates and other complex electronic materials. 



High resolution microscopies of many classes of complex electronic materials such as 
cuprates, manganites, ferroelastic martensites, and relaxor ferroelectric titanates ^[21E1E1 
have revealed previously unsuspected, and puzzling, multiscale modulations of charge, spin, 
polarization and strain variables in stripe- or droplet-like patterning over 1-100 nm scales, or 
up to hundreds of lattice spacings. These sign-varying inhomogeneities or 'texturings' funda- 
mentally affect local electronic, magnetic and structural properties, and appear to be intrinsic: 
arising from the coupling between degrees-of-freedom, rather than induced by extrinsic disor- 
der. 

Ferroelastic martensitic alloys like FePd, and the more complex cuprates and mangan- 
ites, typically have displacive first-order structural transitions PP, described by symmetry- 
adapted strain tensor components as order-parameters (OP), with 1-10 nm criss-cross 'tweed' 
strain variations above, and 10-100 nm 'twins' below, the structural transition pp. There 
is evidence for lattice/charge/spin coupling in both cuprates and manganites in: isotope 
effects doping-dependence of transition temperatures £Q; polaron signatures with 0] 
temperature- and field-dependent four-lobe signatures in diffuse X-ray and neutron scattering; 
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field-dependent wandering ribbons of high conductivity ^1121101 ' anc ^ (fluctuating or static) 
stripe bubbles of modulated buckling-angle and spin [3] . A strain-magnetization coupling is in- 
dicated by cross-responses [7J: e.g., manganites that show colossal magneto- resistance (CMR) 
in magnetic fields also show colossal 'stresso-resistance' (CSR) under hydrostatic pressure. 

The central questions are conceptual: what excitations result, when doped local charges in 
a (para- or antiferro-) magnetic background couple to a nonlinear elastic lattice that supports 
a structural transition? What is the origin of multiscalc intrinsic inhomogcncities? 

We adopt a Ginzburg-Landau (GL) approach, with three generic model assumptions: (1) 
The free energy is harmonic in the non-order-parameter strain components and nonlinear 
in the order-parameter strain (with multiple wells just above, and below, structural tran- 
sition temperatures); (2) The charge density and magnetization variables couple locally, in 
symmetry-allowed ways, to different symmetry-adapted strain-tensor components, (and can 
then act like local internal stresses or temperatures); (3) The charge carriers deforming a unit- 
cell do not damage the lattice by generating defects such as dislocations or vacancies (so there 
is a smoothly compatible matching to the farther, and decreasingly-strained, unit-cells). Long- 
range isotropic Coulomb potentials, with or without, or with electron-phonon coupling, have 
been invoked to explain nanoscale charged stripe patterns 0. However, our central point here 
is that charges and spins locally coupled to strains of a nonlinear lattice, can induce multiscale 
mutual texturings through anisotropic long-range strain-strain forces (themselves Coulombic 
at a chemical-bond level [5]), that arise from generic elastic compatibility constraints |1U) . 

GL model: 1. Strain Free Energy: We consider a two-dimensional (2D) first-order square- 
to-rectangular transformation as a surrogate for 3D tetragonal-orthorhombic PP structural 
transitions; extensions to 3D are possible The Cauchy strain tensor E_ linear in the 

displacement u, is E^ v = (l/2)(A M Uy + A^u^), where A^ is a discrete derivative in the fi = x,y 
directions, and small 'geometric' nonlinearities are ignored. The symmetry-adapted strains are 
the 'rectangular' or deviatoric e = (l/^/2)(E xx — E yy ) , compressional ei = (1 / \/2)(E xx + E yy ) , 
and shear ei = E xy , strains, respectively. The invariant free energy is 

F = (e, ex , e 2 ) + F& (to) + F^ couplm ^ (n, m, e, ei), (1) 

with all quantities scaled |1 lj to be dimensionless. The strain contribution is 

= (a /2) J2A^ £ ) 2 + ^o(e) + F cs (ei,e2), where ^/oq is a strain variation length. The 
Landau term is sixth-order in the deviatoric OP strain Fq = (r — l)e 2 + e 2 (e 2 — l) 2 . Here 
t(T) = (T — T SC )/(T S — T sc ) is a scaled temperature, and for 4/3 > r > 0, Fq has triple wells 
that become degenerate at r(T = T s ) = 1, reducing to double wells for T < T sc or (r < 0). 
The compression/shear (cs) terms are harmonic, F cs = J2p 5^i e ?(r) + k^e^r). 

2. Charge and Magnetization Free Energy: The free energy for a magnetization variable 
m(f) is F^(m) = £ P /[(^)(Am(f)) 2 + (T - T cm )m 2 (r) + |m 4 (r)] - hm. Here / is a 
magnetic/elastic energy ratio, h {^Ja m ) is a magnetic field (length), and the zero-doping 
magnetic transition temperature T cm < T s . The symmetry-allowed couplings are 

F (coup ' m9) (n,m,e,ei) = A n£ ne 2 + A nm nm 2 + [A nl n ei + A ml m 2 ei +p 1 (r)e 1 ], (2) 

r 

where pi(r) > is an external compressional stress. The A n£ ne 2 term is like a local tempera- 
ture ~ re 2 , and for A ne > 0, the charge favours the high-temperature e = unit-cell symmetry. 
The A m im 2 ei term acts like a local stress (~ m 2 ) for e±, and like a local temperature (~ ei) 
for to 2 . We choose A ml > 0, favouring to 2 when spins are closer (ei < 0). The m 2 coeffi- 
cient in F defines an effective local temperature deviation r cm (T, r) = f[T — T c e //(r)]/A ml = 
e i + [.f(T~T cm ) + A nm n]/A m i = Sei(T, r). This implies a temperature-dependent percolation: 



A.R. Bishop et ai: Elasticity-driven Nanoscale Texturing 



3 



on warming, the compressed m 7^ regions below transition ( r cm (T, r) = 5e\{T,r) < 0) will 
shrink. We now mimic manganites and cuprates by further model choices of parameter signs 
and interpretations of variables. 

(i) For 'manganites', with a zero-doping nonferromagnetic T cm = parent compound, 
we focus for simplicity on the magnetization-inducing mobile electrons produced by doping, 
of concentration x e —< n > where the number density n(r) = ^^k 2 / 2vr)e- K l^l is a 
sum over normalized single-charge profiles at sites i, with 2/k the intersite tunneling length. 
Here the A nm nm 2 term with A nrn < means that mobile electrons lock onto and align the 
ferromagnetic m{r) core spins (mimicking a double exchange/ Hund's rule, term). The mobile 
electrons can induce an effective local [B] Curie temperature, T£ff(f) = (— A nm n— A m jei)// > 
0. Since metallicity shrinks unit-cell volume, we take the A nX ne\ term with A n i > 0: n(r) 
favours compression, ei < 0. 

(ii) For 'cuprates', n(f) is the local hole number density of doping fraction Xh =< n > 
into an antiferromagnetic parent compound with Neel temperature T cm ^ 0. With A nm > 
in Anmnm 2 , the staggered magnetization m(r) is reduced by hole doping. As the structural 
transition temperature also decreases rapidly with doping yQ, we take A ne » 1. Since holes 
repel lattice ions, the A n \ne\ term has A n \ < 0: n(f) favours expansion, e\ > 0. 

3. Elastic Compatibility and Anisotropic Long-range Potentials: The St. Venant's com- 
patibility condition JU| expresses the no-defect constraint as A x (A x E_Y = 0. (This 

is analogous to a A.B = no-monopole condition in electromagnetism.) In 2D we have 
|T2| A 2 ei - V8A x A y e 2 = (A 2 X - A 2 y )e, linking OP and non-OP strains. Minimizing the 
free energy F cs with respect to ei,2 while maintaining the compatibility constraint, we find 
e 12 (fc) = B h2 (k)e(k), where S x (fc) = \k 2 - k y 2 ]k 2 /[k 4 + (8A 1 /A 2 ){k x k v ) 2 } and B 2 (k) = 
— (Ai / A 2 )[^/8k x k y / k 2 ]Bi(k) . Substituting back, the seemingly innocuous harmonic F cs yields 
|llllT2*| the crucial ALR potential, encoding unit-cell fourfold symmetries, 

F cs = (A 1 /2)Y / U(kMk)\ 2 ; U(k) = [{k x 2 - k y 2 ) 2 /[k 4 + {%A X /A 2 ){k x k v )% (3) 

k 

while (2) yields nonlocal couplings to the OP. Here U(k) depends on the direction k and 
clearly favors k x = k y diagonal strain textures, with a Meissner-like e-y t2 = expulsion 11 . In 
coordinate space, with f.r' — cos(8 — 8'), the potential U(r — f) ~ cos4(6* — 0')/\f— f 1 ] has 
sign- variation supporting elastic frustration, with the OP strain at a point receiving conflicting 
("ferro/antiferro" ) instructions from other surrounding strains. The power-law (D = 2) decay 
arises from U(k) being scale-free (|fc|-independent) at long wavelengths, rather than from 
proximity to some critical point. 

Our central physical idea is quite simple. Suppose, among (symmetry-broken) rectangular 
unit cells e(f) — 1, that a single unit-cell is made square, e{r) = 0. To maintain lattice 
integrity, the neighboring (and further) unit cells must also deform, with an admixture of non- 
OP strains. As shown in simulations of a nonlinear-strain model under local external stress, 
for Ai/2. >> 1 the large non-OP energy costs can make it profitable to locally summon up the 
(degenerate) competing structure, in energy-lowering higher elastic multipoles: a process of 
adaptive elastic screening |llj . Thus a charge, acting as a local internal stress, can produce an 
unusual sign-varying (i.e. textured) extended polaron modulated by the anisotropic U{r — r'), 
with coupled fields like m(r) also sign-varying. This polaronic elasto- magnetic texture, or 
l pemtoT^ , arises from compatibility and competing ground states, and differs from the more 
familiar magneto-elastic polaron .13; that deforms a single lattice structure. 

We now choose parameters. For the martensite FePd, physical values can be scaled 
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[TT] to be dimensionless, and e — e% = e 2 = 1 correspond to strains ~ 0.02; the scaled 
stress pi = 1 (magnetic field h = 1) corresponds to ~ 0.02 GPa ( ~ 0.25 Tesla); and the 
non-OP and OP elastic constants are A\ = 150, A 2 — 300. We take as illustrative, A\ — 
50, A 2 = 105, ^ao = 0.5, = 1, K = 2, T s = 1, T sc = 0.8, / = 0.3, A ml = +5, and specific 
'manganite' ('cuprate') model parameter sets as T cm — 0, A nm = — l,A n i = +5, A n£ = +2 
{Tcm = 0.6, A nm — +9,^4„i = —5,A ne — 20). This is a regime of globally weak magnetism, 
relatively strong electron-phonon coupling, and dominant compatibility forces. 

Simulation of Textures: The free energy F — F(n, m, e) minimum is found by the 
overdamped limit of a general ferroelastic dynamics |11) : 

• - • - 

de ' dm 

We show selected relaxed profiles ^3] of e, ei,m , with both fc and r plots needed for a full 
understanding. Figure 1 shows strain plots in coordinate space (of ei(r)) and in Fourier space 
(of |e(fc)| 2 ) due to a single charge for 'manganite' parameters. The ne\ local stress term 
would by itself produce a bare single-sign strain, so adaptive elastic screening is responsible 
for the observed textured polaron or 'pemton'. The butterfly- like quadrupolar (and essentially 
cancelling) lobes of both signs in strain and magnetization, extends over ~ 20 lattice spacings, 
explicitly illustrating the concept ^DEIEI of magnetic and structural [2] 'nano-scale phase 
separation'. Bi-pemtons from nearby charges form stripe-like segments. 

For increasing random doping at T — 0.5, the average magnetization < m > rises sharply 
from zero for x e > 0.13, to e.g., < m >= —0.21 at x e — 0.15, with symmetry-breaking 
in the m(r) carried by the pemtons. Fig. 2 shows the mutually deforming multi-pemtons, 
forming meandering ribbons of expanded/compressed unit-cell strain e±(r) or oriented spins 
m(r). Fourier space plots, e.g. of |e(A:)| 2 (or |ei(fc)| 2 ), show a four-lobe shape (as in the 
single-pemton case) reminiscent of diffuse X-ray and neutron scattering |4lll4) . A signature of 
compatibility forces is the k — » — fc inversion symmetry (in the squared strain) relating most 
of the even finescale crinkles. Fig. 2 also shows a reduction, on warming, of both < m > and 
of m(r) percolation |2"lll5|. at fixed doping. 

Figure 3 shows the effects on the composite multipemton of an external spatially varying 
compressional stress pi(r) or magnetic field h(r) , with four quadrants in each picture from the 
long wavelength modulation. The diagonal plots show the p\ — > e\, h — > m or direct responses. 
The off-diagonal plots show the p\ — > m, h e± or cross responses , e.g. fine cloud-like gra- 
dations reminiscent of STM/TEM pictures [2]; and enhanced magnetic percolation under 
stress [3- The CMR/CSR 'colossal' effects can be understood through the "compressed ~ 
magnetic ~ metallic" interconnections, with fields/stresses locally tipping the delicate bal- 
ance between opposite-sign large values (see vertical scales) of m{r),e\(f), thus opening up 
conducting channels both in magnetization (via double exchange), or in compression (via en- 
hanced tunneling). The rich phase diagram J] of charge- (or pemton-) ordered states could 
arise from the orienting long-range compatibility forces ^Hj- Charge profile relaxation into 
locally compressed regions would describe orbital ordering, or coexisting localized/extended 
electronic states. 

For 'cuprate' parameters, the Neel temperature is T cm = 0.6, and at T = 0.5, the parent 
compound has a uniform m. Figure 4 shows that a 'cuprate' single pemton is smaller, stronger, 
and sharper than Fig.l. The 'cuprate' multi-pemton with Xh — 0.1 aligns e(r) into diagonal 
parallel ribbons, forming bubbles of stripes 0, with changes for higher doping. 

In summary, the central insight from our model is that under doping perturbation , a 
nonlinear lattice can produce 'intrinsic inhomogeneities', that are not quenched-defect random 
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spots, but rather, self-organized annealed-texture responses, induced by the multiscale effects 
of local lattice-integrity constraints. These composite textures vary with T,p\ or h. 

The compatible, inter-ce\\ large-strain texturings must be supported by intra-ce\\ deforma- 
tions of the atomic bases ('microstrain'), reflected in bond angle/length distributions |.''>I9| . and 
so will be relevant for complex electronic oxides, with atomic bases of tiltable perovskite octa- 
hedra (that have directionally bonded transition- metal ions, and polarizable/deformable oxy- 
gens .1.) . Further electronic structure studies would allow for both electronic and ionic opti- 
mizations. Further theoretical work includes exploring parameter space extensively, in 2D and 
3D, in overdamped or other JT] regimes; adding charge-hopping dynamics and charge-profile 
relaxations; and strain-related microscopic modelling e.g. of plane bucklings |17| induced 
by octahedral tilts. Further experimental work should include STM mapping of symmetry- 
adapted strains to seek pemton signatures. 

It is a pleasure to thank Seamus Davis, Carlos Frontera, and Venkat Pai for useful conver- 
sations. This work was supported by the U.S. Department of Energy. 
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Fig. 1 - Single polaromc elasto-magnetic texture or pemton state for 'manganite ' parameters: Temper- 
ature T — 0.5. All rows are read left to right. Top row: relief and contour plots of the compressional 
ei(r) strain xlO -3 ; and Fourier space deviatoric strain of |e(fc)| 2 . Bottom row: colour plot of ei(f) 
for separated and nearby pemtons. Note that bi-pemtons form stripe segments. 
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Fig. 2 - Multipemton states for 'manganite' parameters: Electron doping fraction x e = 0.15. Re- 
lief/contour plots; all rows read left to right. Top row: T = 0.5. Plot of the compressional strain 
ei(r). Plot in Fourier space of multiscale |e(fc)| 2 (log scale). Bottom row: magnetization m(f) for 
T = 0.5 and 0.8, with large values compensating to yield < m >= —0.21 and —0.06, respectively. 
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Fig. 4 - Pemtons for 'cuprate' parameters: Temperature T = 0.5; parent compound (subtracted) has 
Neel temperature T cm — 0.6. All rows read left to right. Top row: Colour plot of single-pemton ei(r). 
Relief /contour plot of staggered magnetization m(f) with < m >= 3 x 10~ 3 . Bottom row: Colour 
plots of e(f) for random pemtons with hole fractions x% = 0.1 and 0.2. 
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